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ABSTRACT 

We have decomposed the broad Ha, H/3 and H7 lines of 90 Active Galactic Nuclei 
Q ■ ( AGNs) into a superposition of a very broad and an intermediate Gaussian components 

S— I ' (VBGC and IMGC) and discovered that the two Gaussian components evolve with 

FWHM of the whole emission lines. We suggest that the VBGC and the IMGC are 
produced in different emission regions, namely, Very Broad Line Region (VBLR) and 
Intermediate Line Region (IMLR). The details of the two components of Ha, H/3 and 
H7 lines indicate that the radius obtained from the emission line reverberation map- 
ping normally corresponds to the radius of the VBLR, but the radius obtained from 
t — I the infrared reverberation mapping corresponding to IMLR, i.e., the inner boundary of 

i the dusty torus. The existence of the IMGC may affect the measurement of the black 

hole mass in AGNs. Therefore, the deviation of NLSls from the M-sigma relation may 
be explained naturally in this way. The evolution of the two emission line regions may 
be related to the evolutionary stages of the broad line regions of AGNs from NLSls 
to BLSls. Other evidences for the existence of the IMLR are also presented. 

: 

, Key words: line: profiles; quasars: emission lines; galaxies: structure; Galaxy: evo- 

lution; galaxies: active; galaxies: nuclei 



INTRODUCTION well as the feeding, fueling and growth of supermassive black 

holes. In this paper we attempt to address both problems in 



a coherent way by studying a sample of SDSS AGNs with 
strong broad emission lines. 



Type I AGNs are often classified into two subclasses accord- 
ing to the Full Width at Half Maximum (FWHM) of their 
broad H/3 lines, i.e NLS1 and BLS1. There is a tight cor- 
relation between the black hole mass and the stellar veloc- 
ity dispersion of the bulge of normal galaxies, the so called LINE DECOMPOSITION AND 

M-sigma relation (e.g. Magorrian et al. 1998; Tremaine et ~ STATISTICAL ANALYSIS 
al. 2002;). BLSls are found also to follow this relation well 

(Greene & Ho 2005). However, NLSls seem to deviate from The sample contains 90 objects with clear Ha, H/3, and H7 

this relation; they seem to have much smaller masses or line profiles that can be decomposed. They are selected by 

much higher stellar velocity dispersions (e.g., Wang & Lu La Mura et al (2004) from SDSS 3 based on their balmer 

2001; Zhou et al. 2006). We call this under-massive black line intensities. 21 of them are classified as NLSls. A dou- 

hole problem of NLSls. It is not clear if the luminosity de- ble Gaussian component model, i.e IMGC and VBGC can 

pendence of dust torus and broad line region is simultane- fit the broad line very well. Ha, H/3 and H7 lines are fitted 

ously related to the under-massive black hole problem and with their FWHM fixed with each other; a deviation of 10% 

the receding torus model that used to explain the luminosity are allowed. Several decomposition samples are presented in 

function of AGN (Simpson 2005). Progress towards solving Fig.l. The statistical analysis are presented in Fig. 2. Ha, H/3 

these problems may shed new lights to the understanding and H7 lines behave similarly. Our analysis is focused on H/3. 

the hierarchical evolution of AGN and its host galaxy, as The FWHM ratios range from 1.5 to 4 in the NLS1 group 
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Figure 1. Decomposition example of broad H/3, Ha and H7 lines. 
From (a) to (c), FWHM of the line increases, and the IMGC 
becomes weaker and may disappear sometimes, e.g., in subpancls 
cl & c3. Note that for in subpanel c2, the IMGC of Ha line 
is still detectable, in contrast to H/3 and H7 lines, because the 
intermediate component of Ho is normally much stronger than 
H/3 and H7 lines. 



which lie on the left of Fig.2(al). The intensity ratio of the 
very broad Gaussian component to the whole line is about 
0.6 in these objects. This means that for NLSls the inter- 
mediate Gaussian component has an intensity comparable 
to the very broad one, and consequently the FWHM of the 
whole line of NLSls is dominated by the intermediate Gaus- 
sian component. With FWHM increasing, the intermediate 
Gaussian component becomes weaker and finally disappears, 
which can be seen from Fig.2(bl). The whole line can be 
simply described by a single very broad Gaussian compo- 
nent when FWHM reaches about 5000 km-s" 1 . In Fig.2(cl), 
FWHM ratio of the intermediate Gaussian component to 
the very broad Gaussian component becomes larger with 
increasing FWHM. 
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Figure 2. Statistical analysis of broad H/3, Ha and H7 lines. 
With the FWHM increasing, FWHM ratio of VBGC to the whole 
line becomes smaller and finally reaches unity, the intensity ratio 
of the VBGC to the whole line becomes larger and finally also 
reaches unity, and the IMGC becomes broader and weaker. 
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3 BHM CORRECTION 

Reverberation-based black hole mass is calculated by the 
virial equation 

Rblr. 



Mbh = 



G 



:/ 2 FWHM 2 



lid . 



(1) 



where FWHMh/3 is the FWHM of the whole H/3, which rep- 
resents the virial velocity of the BLR. For an isotropic veloc- 
ity distribution, as generally assumed, / = y/3/2. Rblr is 
the BLR radius obtained by reverberation mapping, which is 
related to an AGN's continuum luminosity by the empirical 
Radius-Luminosity relation (Bentz et al. 2006): 



los 



-Rblr 
It — days 



(0.518±0.039) lo{ 



AL A (5100l) 



-21.2±1.7.(2) 



erg • s 

This relation has been starlight-corrected. We use this equa- 
tion to calculate Rblr throughout this work. Because re- 
verberation mapping is used to measure the radius of the 
broad line region, naturally it normally measures the radius 



Figure 3. Black hole masses (in units of solar mass Mq) are 
plotted in logarithm scale. Circles are black hole mass calculated 
with equation (1). Dots are that after correction using equation 
(3). Clearly the correction is more effective for AGNs with smaller 
FWHM, i.e., NLSls. 



of the inner most emission line region, i.e. the VBLR in our 
analysis. Consequently, FWHM of the very broad Gaussian 
component, instead of FWHM of the whole line, should be 
used to calculate the black hole mass. We therefore correct 
the black hole mass in this way: 

«»^ f ™M-,=f, ,3, 

where FWHMb is the FWHM of the very broad Gaussian 
component, and Rvblr is taken as the 7?blr. in the R-L re- 
lation. It gives a more significant mass correction for NLSls 
than for BLSls as shown in Fig. 3. After such correction, 
NLSls still have smaller black hole masses but normal ac- 
cretion rate in units of the Eddington rate. Fifteen objects 
in our samples have velocity dispersion measurement data 
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(sigma) in Shen et al (2008). The masses of all NSL1 are well 
below, but become more very close to, the predictions of the 
M-sigma relation before and after the correction. The effect 
of the correction is obvious, albeit small number statistics 
due to the limited sample. 



4 THE LOCATION OF THE IMLR AND THE 
EVIDENCES FOR ITS EXISTENCE 

4.1 Broad Line Region evolution 

We assume that both of the two regions (VBLR and IMLR) 
are bounded by the central black hole's gravity. Then, we 
have 



ft 



RvBLR 



GMbh 
R 2 

VBLR. 
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These two equations lead to 



RlMLR. = ( 



fl VvBLR. 
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where Co = (jj-) 2 . Co is a constant which need to be deter- 
mined by further studies. The radius measured from rever- 
beration is taken as Rvblr- Therefore, we can calculate the 
radius of IMLR from the above equation. Rimlr is propor- 
tional to Co- Co = 1 is used for all the calculation below. 
fiiMLH is obtained from Ha and H/3 lines separately. Fig. 4 
is the evolution of -Rvblr and -Rimlr with luminosity. The 
radius of the IMLR increases slower than VBLR with lu- 
minosity and black hole mass increasing. The two emission 
regions have a trend to merge into one with larger higher 
luminosity and black hole mass, as -Rimlr oc L5100°' 37±o ' 06 
while Rvblr oc L5100 a52±0 ' 04 . 



4.2 The location of IMLR 

The relation Rimlr oc L5100°' 37±o ' 06 is consistent with the 
receding velocity of torus based on infrared reverberation 
mapping (Suganuma et al. 2006) and the receding velocity of 
torus that we can calculate based on analysis of type I AGN 
fraction (Simpson 2005). Our other analysis also shows that 
the IMLR may have higher density, more dusty than the 
VBGC, and the IMGC show slightly Baldwin effect while 
the VBGC does not. The Baldwin effect can be explained 
as a flattened geometry. These results suggest that IMLR 
is the inner hot skin of torus. A simple scenario of the 
BLR (BLR=VBLR+IMLR) evolution can be constructed 
as shown in Fig. 5. The inner spherical region is the VBLR 
(Very Broad Line Region). It expands to a larger radius 
with luminosity increase and perhaps also black hole mass 
increase. IMLR is the inner part of the torus, which can 
be sublimated by the central radiation and thus its radius 
also increases with luminosity increase. When the luminos- 
ity is high enough, the broad line region and torus become 
one single entity, but with different physical conditions. This 
scenario is consistent with the dust bound BLR hypothesis 
which also showed that the delay time of infrared emission 
is always longer than the delay time of emission lines of 
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Figure 4. Correlation of the radii of IMLR and VBLR with lu- 
minosity. Circles are the VBLR radius, X represent IMLR, radius 
derived from H/3 lines and + from Ha lines. Circles filled with dots 
represent the objects missing the intermediate Gaussian compo- 
nent. The radius of the VBLR is obtained from the i?BLR. ~L5100 
(equation 2 not equation 3). These correlations supports a sce- 
nario of hierarchical evolution of AGNs. 




Figure 5. Cartoon of the Broad Line Region evolution. Contin- 
uum photons from the central region have slightly higher energies 
along the nearly edge-on direction. With black hole mass and lu- 
minosity increasing, both the VBLR and IMLR expand. The ra- 
dius of VBLR increases faster, so the two regions have a trend to 
merge into one. 



the corresponding AGN, but sometimes they are almost the 
same (Laor 2004; Elitzur 2006). Suganuma et al. (2006). 



4.3 Evidence for the existence of IMLR 

The first concrete evidence comes from the the study of the 
micro-lensing of the Broad Line Region in the lensed quasar 
J1131-1231 (Sluse et al. 2007, 2008). In this work they found 
evidence that the H/3 emission line (as well as Ha) is differ- 
entially microlensed, with the broadest component (FWHM 
4000 km/s) being much more micro- lensed than the nar- 
rower component (FWHM 2000 km/s). Because the ampli- 
tude of micro-lensing depends on the size of the emitting 
region, this is a clear evidence that the broadest compo- 
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Figure 6. Different responses of different geometries to a delta 
impulse. The VBLR is assumed to be a spherical shell of uniform 
density, and the IMLR is assumed to be a cylindrical shell of 
uniform density. Time is normalized to arbitrary unit. Each figure 
shows a different -Rimlr — ^?VBLR pair. Clearly in all cases the 
responses of IMLR have shorter variability and thus will result in 
narrower RMS spectra than that from VBLR. 



nent of the emission line is more compact than the narrower 
component. 

Then, we turn to the reverberation mapping experi- 
ment. If the IMLR is really a discrete emission region apart 
from the VBLR, another peak corresponding to the radius 
of this region might appear in the cross correlation function 
between the lightcurves of the emission line and the con- 
tinuum. In fact, a possible example exists in the database 
of reverberation mapping observations (Peterson et al. 1998; 
2004) . Mrk 79 shows obviously double peaks in the cross cor- 
relation centroid distribution which has not been explained 
well so far which can be well explained in our model (please 
refer to Zhu et al. 2008 for details). 

Another fact known to us is that the RMS spectrum of 
the broad emission line is normally (but not always) nar- 
rower systematically than the mean spectrum (Collin et al. 
2006). If the emission lines are produced from two distinct 
regions, VBLR and IMLR, they will respond to the central 
continuum radiation differently. Because the IMLR has a 
flattened geometry, it can create a narrower response func- 
tion to a delta impulse even with radius much larger than 
the VBLR, which is confirmed by our calculations shown in 
Fig.6. 



5 SUMMARY 

We conclude that the broad LL3, Ha, and H7 line of the 
AGNs in our sample can be well fitted by a double Gaus- 
sian model, i.e., VBGC and IMGC. The properties of the 
two components indicate that the two Gaussian components 
originate from two emission regions, i.e. VBLR and IMLR. 
The physical properties of IMLR suggest it should be the 
torus inner boundary region. Based on our model, The black 
hole mass measured by RM should be corrected, after cor- 



rection the NLS1 follows M-sigma relation as well as BLS1. 
Therefore the problem of under-massive black hole in NLSls 
appears to have disappeared. Our results offer strong evi- 
dence of the evolution of the broad line region (consists of a 
IMLR and a VBLR) from NLSls to BLSls in the hierarchi- 
cal evolution scenario. 

Finally, we show that the study of micro-lensing provide 
concrete evidence for the existence of IMLR and Mrk 79 pro- 
vides possible evidence. Narrower RMS spectrum also sup- 
port our model. The existence of a weak intermediate Gaus- 
sian component in the broad emission lines of many sources 
with reverberation mapping measurements may cause sys- 
tematic biases for the measurements of the radius of the 
VBLR, It will be helpful to decompose each broad emission 
line into two Gaussian components as we have done here, 
and then do cross-correlation analysis between the contin- 
uum and each of the two components, in order to measure 
the Radius-Luminosity relations for the two components in- 
dependently. This would provide decisive test for our model 
presented in this paper, as well as providing more accurate 
measurements of the masses of supermassive black holes in 
AGNs. 
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